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We study an anomalous behavior of the height fluctuation width in the crossover from random 
to coherent growths of surface for a stochastic model. In the model, random numbers are assigned 
on perimeter sites of surface, representing pinning strengths of disordered media. At each time, 
surface is advanced at the site having minimum pinning strength in a random subset of system 
rather than having global minimum. The subset is composed of a randomly selected site and its 
(£— 1) neighbors. The height fluctuation width W 2 (L;£) exhibits the non-monotonic behavior with 
I and it has a minimum at £*. It is found numerically that £* scales as t ~ L 0,59 , and the height 
fluctuation width at that minimum, W 2 (L;£*), scales as ~ L ' 85 in 1+1 dimensions. It is found 
that the subset-size I* (L) is the characteristic size of the crossover from the random surface growth 
in the KPZ universality, to the coherent surface growth in the directed percolation universality. 
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Recently subjects in the field of non-equilibrium sur- 
face growth have been attractive, which is due to the 
interest in theoretical classification of universality and 
also due to their applications to various physical phe- 
nomena such as crystal growth, molecular beam epitaxy, 
vapor deposition, and biological evolution, etc pj. The 
interesting feature of non-equilibrium surface growth is 
the nontrivial scaling behavior of the height fluctuation 
width j2j, i.e., 

W 2 (L, t) = (A - hf) ~ L 2a f{t/L% (1) 

i 

where hi is the height of site i on substrate. Here, 
h, L, and d! denote mean height, system size, and sub- 
strate dimension, respectively. The symbol (• • •) stands 
for statistical average. The scaling function behaves as 
f(x) — > constant for x >■ 1, and f(x) ~ x 2lB for iCl 
with (3 = a/z. The exponents, a, (3 and z are called 
the roughness exponent, the growth exponent, and the 
dynamic exponent, respectively. 

The Kardar-Parisi-Zhang (KPZ) equation || was in- 
troduced to account for the effect of sideway growth, 
which is written as 



dh(x, t) 
dt 



vV 2 h{x,t) + ^(V/i(x,i)) 2 + 77(35, t), (2) 



where rj(x,t), called thermal noise, is assumed as white 
noise, (r)(x,t)) = and (77(2;, t)r)(x', if)) = 2DS d ' (x - 
x')6(t — t') with noise strength D, Many stochastic mod- 
els in the KPZ universality class have been introduced 
[jij. Among them, the restricted solid-on-solid (RSOS) 
model was introduced by Kim and Kosterlitz jj] , satisfy- 
ing the scaling relation, at + z t = 2. The subscript means 
that the exponents are for thermal noise. In the RSOS 
model, a particle is deposited at a randomly selected 
site as long as the height difference Ah between nearest 



neighbor columns remains as Ah < 1 even after deposi- 
tion, otherwise, the particle is excluded from deposition. 
On the other hand, one may modify the dynamic rule of 
the RSOS model by replacing the exclusion by avalanche: 
At each time step, a particle is deposited at a randomly 
selected site, and avalanche may occur successively to 
nearest neighbor sites as long as the height difference 
between nearest neighbors, Ah > 1. Then it was argued 
that this model also belongs to the KPZ universality class 
Q , but the model requires relatively large system size to 
see its asymptotic behavior. Let us call the former model 
the RSOS-A model, and the latter model the RSOS-B 
model. The surface-growth in the KPZ universality class 
is called the random surface-growth. 

Physical properties of growing surface in disordered 
media are different from those of the thermal KPZ equa- 
tion, Eq. (2). In order to account for the effect of disorder 
in porous media, quenched noise which depends on 
position x and height h, replaces thermal noise in Eq. (2). 
Then the quenched KPZ (QKPZ) equation is written as 



dh(x, t) 
dt 



v\/ 2 h(x, t) + ^(Vft(x, t)) 2 + r)(x, h), (3) 



where the noise satisfies, {rj(x.h)) = 0, and 
(r](x,h)ri(x',ti)) = 2D5 d ' (x - x')5(h - h'). Stochastic 
models associated with the QKPZ equation have been 
introduced J7U8J. The models show that the surface of 
the QKPZ equation in 1+1 dimensions belongs to the 
directed percolation (DP) universality. The roughness 
exponent a q in the QKPZ equation is given by the ratio 
of the correlation length exponents in perpendicular and 



parallel directions, v± and v 
which is a n 



II , of directed percolating 
0.63. On the other 



clusters, which is a q = v±/v\\ 
hand, recently Sneppen introduced a stochastic model 
in which surface grows coherently ||. In that model, 
random numbers, representing disorder of porous media, 



are assigned on each perimeter site of surface. Surface 
is advanced at the site having global minimum among 
the random numbers. The avalanche rule is then applied 
successively to nearest neighbor sites as long as A/i > 1. 
Random numbers at the columns with increased heights 
are updated by new ones. The Sneppen model also be- 
longs to the DP universality, and the roughness exponent 
is a s « 0.63 in 1+1 dimensions The surface-growth of 
the Sneppen model is called the coherent surface-growth. 

The coherent surface-growth model is closely related 
to the self-organized evolution model for biological sys- 
tems [jLOj . In the evolution model, one considers random 
numbers assigned in one-dimensional array, which rep- 
resent fitness of each species. Mutation of species is 
described in the model by updating random numbers. 
The updating occurs at the site having global minimum 
random number and its two nearest neighbors, and the 
random numbers at those sites are replaced by new ones 
at each time. Then as times go on, relatively small ran- 
dom numbers disappear by the updating process, and the 
distribution of random numbers exhibits a self-organized 
critical behavior. On the other hand, one may think of 
the situation where biological evolution is not motivated 
in globally optimized manner, but it may be driven by 
optimization within finite region out of entire system. 
This is analogous to the case that spin glass system is 
in meta-stable state within finite relaxation time rather 
than in globally stable state. Motivated by this idea, in 
this paper, we introduce a surface growth model in dis- 
ordered media, where surface growth occurs at the site 
having the minimum random number in a subset of the 
entire system rather than having the global minimum 
random number. We think this model might be relevant 
to the case where the relaxation of surface growth in dis- 
ordered media is not fast enough to spread into the whole 
system, so that surface growth is driven not in globally 
optimized manner, but in locally optimized manner. 

To be specific, the model we consider in this paper is 
defined as follows: First, we consider one dimensional 
flat substrate with system size L. Random numbers are 
assigned on each site, which represent energy barriers 
(fitness) in the evolution model or pinning forces in the 
surface growth model by Sneppen. Second, we select 
a site randomly, and consider the subset composed of 
£ = 2r + 1 elements, the randomly selected site and its 
2r neighbor sites within distance r. It is worthwhile to 
note that the subset is regarded as a random sample 
because the site in the middle of the subset was selected 
at random. The subset is formed instantaneously, and 
its territory might overlap with subsequent one as shown 
in Fig. 1. Next, surface is advanced at the site having 
minimum random number among the £ elements in the 
subset, and the avalanche process is followed successively 
at its neighbor sites to keep the RSOS condition, Ah < 1, 
and it may spread out over the boundary of the subset. 
Finally, the random numbers at the sites with increased 



height are updated with new ones. The dynamic rule of 
the model is depicted in Fig. 1. When 1=1, this model 
corresponds to the RSOS-B model in the KPZ univer- 
sality, whereas when I = L, it does to the Sneppen-B 
model in the DP universality. Accordingly one may see 
the crossover behavior from the KPZ limit to the DP 
limit as increasing I. 

Since the roughness exponent at = 1/2 in the KPZ 
limit is smaller than the one a q « 0.63 in the DP limit, 
one may expect at a glance that the height fluctuation 
width W 2 (L;£) in steady state increases monotonically 
with increasing £. However, we found the anomalous 
behavior numerically that W 2 (L;£) decreases with in- 
creasing £ for small £, and increases for large £ as shown 
in Fig. 2. The minimum of W 2 (L;£) becomes steeper 
and its location, £* /L, which was rescaled by system size 
L, approaches to zero as L increases. It is found that 
the location of the minimum scales as £* ~ L 59 , and 
the height fluctuation width at this minimum scales as 
W 2 (L;£*) ~ L - 85 , which are shown in Figs. 3 and 4. 
The estimated values of £* and W 2 (L;£*) for different 
system sizes are tabulated below. 

Table 1. Numerical values of the location £* and the 
height fluctuation width W 2 (L;£*) at the minimum for 
different system sizes. 



L 


£* 


W' Z {L;£*) 


64 


11 


5.9 


128 


19 


11.2 


256 


27 


20.3 


512 


39 


35.8 


1024 


59 


62.5 


2048 


95 


110.0 



The anomalous behavior may be attributed to the two 
effects, the random effect for small £ and the coherent 
effect for large £. For small £, surface grows by random 
deposition with avalanches, and it belongs to the KPZ 
universality. Thus the height fluctuation width depends 
on system size as W 2 (L;£) ~ L 2at with 2at — 1, how- 
ever, it would also depend on subset size I. In order to 
find out the ^-dependent behavior of W 2 (L;£) in phe- 
nomenological level, we plot W 2 (L;£) versus £ in double 
logarithmic scales for several values of L as shown in 
Fig. 5. The slopes are measured to be w —0.2 for large 
system sizes L — 1024 and 2048. Based on this measure- 
ment, W 2 (L; £) is written as W 2 (L\ £) ~ £~°- 2 L for small 
£. This result is contrary to the one, based on the coarse- 
graining scaling argument, W 2 (L;£) ~ (L / £) 2at £ 2a " , 
which exhibits the increasing behavior of W 2 with in- 
creasing £. On the other hand, when £ is large enough, 
surface updating is initiated mainly at the site having 
global minimum random number of entire system. When 
the site of global minimum is selected, which occurs with 
probability i/L, the surface becomes correlated by the 
Sneppen dynamics, however, when the site of the global 
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minimum is not selected, the correlation formed by the 
Sneppen's dynamics is relaxed. Thus, when £ is large 
enough so that the contribution by the Sneppen dynam- 
ics is sufficiently dominant, one may write the domi- 
nant term of the averaged height fluctuation width as 
W 2 (L;£) — (^L 2 "". That is because the statistical av- 
erage was taken over the quantity of the square of mean- 
height deviation in Eq. (1). Combining the two asymp- 
totic behaviors obtained in phenomenological level, the 
height fluctuation width is written as 

W 2 (L;£)^r°- 2 L+d)L 1 - 26 . (4) 

The two terms exhibit the competing behavior with re- 
spect to £, which yields the anomalous behavior. Taking 
the derivative with respect to £, the location £* of the 
minimum is obtained as that £* ~ L - 62 , which is close to 
the numerical measurement, t* ~ L 59 . We also exam- 
ine the size-dependent behavior of the height fluctuation 
width at the minimum by plugging £* into Eq. (4), and 
obtain that W 2 (L;£*) ~ i 88 . This result is also close 
to the numerical measurement, W 2 (L;£*) <~ L°- S5 . The 
numerical estimations for £* and W 2 (L;£*) are better 
explained by minimizing the formula, 

t>2a t f 

W 2 (L;£)^(L^)L^ + (i)L^, (5) 
however, the derivation of this formula is not clear. 

In order to understand physical meaning of the char- 
acteristic subset-size £*, we plot W 2 (L;£) versus L up to 
L = 2048 in double logarithmic scales for typical subset 
sizes, £ = 5 and 41 in Fig. 6. The size of £ — 5 corre- 
sponds to the case where it is smaller than £*(L) for all 
system sizes L used in Fig. 6. However, the size of £ = 41 
does to the case where it is smaller than £*(L) in part 
for L = 1024 and 2048, close to t{L) for L = 512, but 
larger than £*(L) in part for L = 64 and 128. For the 
case of I = 5, all data are on a straight line, whereas, for 
the case of I = 41, forming a straight line breaks down 
for relatively smaller system sizes L — 64 and 128. Fig. 6 
suggests that £* (L) be the characteristic subset-size such 
that when t < £*(L), the roughness of overall surface is 
determined by random effect, whereas when £ > £*(L), 
it is done by coherent effect. Accordingly, the subset size 
£* has the meanings of not only the location of the mini- 
mum of the anomalous height fluctuation width, but also 
the critical size at which the crossover from random to 
coherent surface-growths occurs. Also note that even for 
the case of £ > £*(L), the surface height fluctuation width 
does not behave as W 2 (L;£) ~ L 2a " . The roughness ex- 
ponent has smaller value than 2a q as appeared in Fig. 6, 
and the numerical estimation of the roughness exponent 
is likely to be 2a q — 1 as Eq. (4) for fixed £. The behav- 
ior of W 2 (L;£) ~ L 2oLq occurs when £ also increases as 
L increases. In thermodynamic limit, L — » oo, the char- 
acteristic size £* goes to infinity, so that for finite £, the 



roughness of overall surface is determined by random ef- 
fect, and the surface is described in the KPZ universality. 
Next, we examine the height- height correlation function, 

C{r-L,£)^(^( h ^ + r)-h(x)) 2 ), (6) 

X 

which is defined for fixed L, and I. In Fig. 7, we plot 
C{r)/l 2aq versus (r/£) in double logarithmic scales for 
several values of £, and L = 2048. In Fig. 7, the data are 
well collapsed for rjl < 1, while they are not collapsed 
for r/£ > 1. Accordingly, the coherent surface growth 
occurs within the range of r < £, however, the roughness 
of overall surface is determined by the criterion depend- 
ing on £* (L) and L above. 

It would be interesting to study the ^-dependent behav- 
ior of dynamic properties of the height fluctuation width 
W 2 (L, t; £). The study is based on numerical simulations 
for fixed system size, say, L = 1024. As shown in inset 
of Fig. 8, there exist four distinct regimes for W 2 (L, t; £). 
In the first regime, W 2 (L,t;£) increases according to 
the Poisson distribution, and W 2 (L,t;£) ~ t 2 ^ 1 with 
2/?i = 1. The first regime terminates at t\, which is 
independent of subset-size £. In the second regime, sur- 
face becomes correlated by the coherent effect, which is 
caused by the selection of minimum random number in 
selected subset, however, the decorrelation also occurs 
simultaneously by the random effect, which is caused by 
the selection of the random subset. Since the value of 
the dynamic exponent z s = 0.63 for the Sneppen dy- 
namics is smaller than the one Zt — 1.5 for the KPZ 
dynamics, the coherent effect spreads faster than the 
random effect in early times. Thus the growth exponent 
2/32 in the second regime has the value more likely close 
to the Sneppen value, 2(3 S = 2, however, the value is 
a little bit smaller by the decorrelation by the random 
effect. The growth exponent 2/3 2 depends on subset-size 
£ as tabulated in Table 2. Based on the measurement in 
Fig. 9, the growth exponent is likely to depend on £ as 
2/3 2 ~ (0.1) log£ for £ < £*, however for £ > £*, the value 
of the growth exponent is expected to be close to the one 
of the Sneppen dynamics as shown in Fig. 8. The second 
regime terminates at t 2 - In Fig. 10, the threshold times 
are likely to scale as t 2 ~ £ ' 37 for £ < £* . The values 
of the height fluctuation width at the threshold value 
t 2 are likely to scale as W 2 {L,t2\£) ~ £ ' 92 as shown in 
Fig. 11. In the third regime, the random effect appears 
much dominantly, and the coherence of surface formed in 
the second regime becomes decorrelated in this regime. 
As subset size is smaller, the third regime is much dom- 
inant, and the growth exponent (3$ becomes much closer 
to the KPZ value, whereas as subset size is larger, the 
decorrelation effect becomes much weaker, so that the 
growth exponent (3% becomes smaller. It is likely that 
2/3 3 ~ (—0.37) log£ as shown in Fig. 12. The third regime 
terminates at t% . The numerical values of t% for different 
sizes of £ locate too closely for small £ to be measured 
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numerically. 

Table 2. Numerical estimation for the values of the 
growth exponents, the threshold times, and the height 
fluctuation widths for various subset sizes I. 



£ 


2/3 2 


!ogio(*2) 


logio(W?) 


2/3 3 


3 


1.12 


-0.10 


0.20 


0.66 


5 


1.16 


-0.10 


0.20 


0.64 


9 


1.23 


0.00 


0.33 


0.60 


17 


1.31 


0.10 


0.55 


0.50 


33 


1.40 


0.20 


0.80 


0.38 


65 


1.59 


0.31 


1.12 


0.28 


129 


1.68 


0.42 


1.60 


0.17 


257 


1.78 


0.60 


1.80 


0.17 



We also investigated the ^-dependent behavior of the 
distribution of random numbers after reaching saturated 
state. As shown in Fig. 13, the distribution is flat for 
£ — 1, and exhibits a critical behavior for £ = L. Between 
the two limits, the distributions look like a rounded step 
function. It would be interesting to note that all distri- 
bution functions for different i pass through a specific 
value of random number B c , which corresponds to the 
threshold of the self-organized critical state M. The 
value of B c equals to 1 — P c = 0.462, where P c is the 
directed percolation threshold. 

Recently, the crossover behavior from the random 
surface growth to the coherent surface growth was 
considered by Vergeles |ll[. In that study, surface 
growth occurs at site x on substrate with the proba- 
bility, P{x) ~ e~~ q ( x ^ T , where q(x) is a random pinning 
strength of site x, and is also updated with height ad- 
vance. T is temperature. It was found that the surface 
of the model reduces to the one of the Sneppen dynamics 
when T = 0, however, for T ^ 0, it does to the one of 
the KPZ universality. The temperature T plays a role of 
tuning parameter for the crossover behavior. However, 
since the tuning parameter is in the form of exponential 
function, it is hard to see finite size dependent behavior 
of the crossover which is very sensitive to tuning the 
parameter as we studied in this paper. Nevertheless, 
the anomalous behavior may be observed barely in the 
plot of W(L; T) versus L for different temperatures, Fig. 
1(b) in Ref. 11, where the curves of W(L,T) cross to 
each other, However the crossover behavior has not been 
remarked in Ref. 11. 

In summary, we have introduced a stochastic model 
for surface growth, which is a generalization of the re- 
stricted solid-on-solid model in the KPZ universality and 
the Sneppen model in the directed percolation limit, and 
have investigated the crossover of the two limits. De- 
position occurs at the site having minimum of random 
numbers within finite subset rather than of entire sys- 
tem. The subset is composed of £ elements, a randomly 



selected site and its € — 1 neighbors. Changing subset-size 
£, the height fluctuation width exhibits the anomalous 
behavior having a minimum. The anomalous behavior 
is due to the two competing effects, the random effect 
for small £ and the coherent effect for large £. The min- 
imum of the surface height fluctuation width, locating 
at I* ~ L - 59 , is scaled as W 2 (L;t) ~ L°- S5 in 1+1 
dimensions. The characteristic subset size £*(L) has the 
meaning that for £ < £*{L), the surface grows randomly 
and belongs to the KPZ universality, whereas for the 
opposite case, surface grows coherently. The dynamic 
properties of the crossover have also been investigated. 
In early stage of growth, surface becomes correlated ac- 
cording to the Sneppen dynamics, and in late stage, the 
surface correlation is relaxed by the random process. The 
phenomenon of the dynamic correlation-decorrelation be- 
havior also appears in a stochastic model Jl2| for the flux 
line dynamics with transversal and longitudinal fluctua- 
tions, which might be described by the coupled quenched 
KPZ equation (u|. 
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FIG. 2. Plot of the surface height fluctuation widths 
W 2 (L;£) versus subset size l/L rescaled by system size 
L in steady state. Numerical data are for system sizes 
L = 64, 128, 256, 512, 1024 and 2048 from the bottom. 




FIG. 1. Schematic representation of the stochastic rule. FIG. 3. Double-logarithmic plot of the estimated location 

The sites arrowed are randomly selected sites. The dark I* of the minimum versus system size L. The data are for 

squares denote the sites of minimum random number within L = 64, 128, 256, 512, 1024, and 2048. 
subset sized i = 5. The white squares denote the sites up- 
dated by avalanches. The subset could overlap with the one 
(lower one) of subsequent time. 
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FIG. 4. Double-logarithmic plot of the surface height fluc- 
tuation width at the minimum position W 2 (L;£*) versus sys- 
tem size L. The data are for L = 64, 128, 256, 512, 1024, and 
2048. 
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FIG. 7. Double-logarithmic plot of C(r)/£ 2aq versus r/£ 
for subset sizes i = 5, 11, 21, 41, 61, 81 and 141. The data are 
for system size L = 2048. 



FIG. 5. Double-logarithmic plot of W 2 {L;£) versus I for 
various system sizes. The data seem to be on straight lines 
with slope —0.2. 





FIG. 6. Double-logarithmic plot of W 2 (L;£) versus L for 
typical subset sizes £ = 5 and 41. 
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FIG. 8. Double-logarithmic plot of W 2 {L, t; £) versus time 
t for system size L = 1024. The data are for subset sizes 
£ = 3, 5,9, 17, 33,65, 129, 257 and 513. Inset: Distinct four 
regimes are observed for the case of L = 1024 and I — 65. 
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FIG. 9. Semi-logarithmic plot of 2/3 2 versus £ for FIG. 11. Double-logarithmic plot of W 2 2 versus I for 

L = 1024. The data seem to be on a straight line with slope L = 1024. The data seem to be on a straight line with slope 
0.1 up to the characteristic subset size £*. 0.92 up to the characteristic subset size £*. 
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FIG. 10. Double-logarithmic plot of ti versus £ for 
L — 1024. The data seem to be on a straight line with slope 
0.37. 
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FIG. 12. Semi-logarithmic plot of 2/33 versus £ for 
L = 1024. The line has the slope —0.37, and the data are 
for I = 3, 5, 9, 17, 33, 65 and 129, respectively. 
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FIG. 13. Plot of the random number distribution after 
saturation for L — 256. The data are for i = 1,3, 5, 27, 129 
and 256, respectively. 



